A path to poor coherence in heavy fermions from Mott physics and hybridization 
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We investigate the anomalous metal arising by hole doping the Mott insulating state of the periodic 
Anderson model. Using Dynamical Mean- Field Theory we show that, as opposed to the electron- 
doped case, in the hole-doped regime the hybridization between localized and delocalized orbitals 
leads to the formation of composite quasi-particles reminiscent of the Zhang-Rice singlets. We 
compute the coherence temperature of this state, showing its extremely small value at low doping. 
As a consequence the weakly-doped Mott state deviates from the predictions of Fermi-liquid theory 
already at small temperatures. The onset of the Zhang-Rice state and of the consequent poor 
coherence is due to the electronic structure in which both localized and itinerant carriers have to 
be involved in the formation of the conduction states and to the proximity to the Mott state. By 
investigating the magnetic properties of this state, we discuss the relation between the anomalous 
metallic properties and the behavior of the magnetic degrees of freedom. 

PACS numbers: 



I. INTRODUCTION 

The rise of the field of strongly correlated materials 
revealed a number of unexpected intriguing phenomena 
which can not be explained within the standard theory 
of solids.^ The paradigm of correlation effects is based on 
the Mott insulating state and the Mott-Hubbard metal- 
insulator transition, ^''^ but a key role is also played by 
high-temperature superconductivity in copper oxides'*'^ 
and the unconventional superconductivity at the edge of 
a magnetic phase observed in heavy fermions. More re- 
cently, the partnership between exotic superconductivity, 
strong correlations and magnetism has been strengthened 
by the discoveries in the iron-based superconductors,^ 
in the alkali-doped fuUerides^'^'^ and possibly also in the 
molecular conductors based on aromatic molecules. ^^"^'^ 

A common companion of Mott physics and anoma- 
lous superconductivity is the deviation from the standard 
Fermi-Liquid (FL) theory in the metallic phase, or non- 
Fermi- liquid (NFL) behavior.^** The FL theory describes 
a system of interacting fermions as a collection of renor- 
malized non-interacting quasi-particles which propagate 
coherently in the solid. The main qualitative effect of 
the electron-electron correlations is to enhance the effec- 
tive mass and accordingly reduce the coherence of the 
Fermi gas. This reflects in a reduction of the coherence 
temperature, the scale at which the thermal fluctuations 
destroy the coherent motion. However in many com- 
pounds, most notably heavy fermion materials and un- 
derdoped cuprates, this picture breaks down and the car- 
riers can no longer be described as long-lived excitations 
as they acquire a finite lifetime. This behavior directly 
influences the transport properties leading to anomalies 
in the temperature dependence of the resistivity. 

In this paper we present a general mechanism based 



on Mott physics and multi-band effects which leads to 
a metallic state with an extremely small FL coherence 
temperature. Empirically, this system will display a NFL 
behavior already at exceedingly small temperatures. The 
key element is the hybridization between a strongly cor- 
related band and a weakly interacting band that leads 
to the formation of hybrid entities. The binding with 
the localized /-electrons hinders the motion of the car- 
riers leading to a coherence temperature which is much 
smaller than the (already renormalized) scale predicted 
by FL theory on the basis of mass renormalization. 

Our approach is based on the periodic Anderson model 
(PAM) , a widely accepted correlated electrons model for 
the description of the heavy fermion physics. In its min- 
imal form the PAM describes a set of non-dispersive 
strongly correlated electrons, hybridizing with a band of 
conduction electrons. In a general framework the PAM 
provides a more detailed description of the electronic 
configuration of correlated materials with respect to the 
Hubbard model, by taking into account the effects of the 
inclusion of non-correlated bands. We solved the PAM 
using dynamical mean- field theory (DMFT),-"^^ one of the 
most powerful and reliable tools to study correlated ma- 
terials. 

Following previous studies^^"^^ we investigate the 
model around the Mott insulating state which takes place 
for large interactions and odd integer total occupation. 
The doping-driven transition has been thoroughly inves- 
tigated in Ref. 17,19, and a NFL behavior in the hole- 
doped side has been demonstrated in Ref. 18. Here we 
extend this work by analyzing the coherence-incoherence 
crossover which leads to the NFL behavior and its de- 
pendence on doping. We will therefore focus on the scat- 
tering properties of the system and we will detail their 
relation with the magnetic degrees of freedom. Finally 
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we establish a connection between the finite-temperature 
breakdown of the FL and the competition between anti- 
ferromagnetic and ferromagnetic short-ranged correla- 
tions. 

The manuscript is structured as follows. In Sec. II we 
introduce the PAM and the related DMFT equations. In 
section Sec. Ill wc briefly discuss the doping-driven Mott 
transition in the PAM. In Sec. IV we present the main re- 
sults of this work, namely the strongly incoherent nature 
of the low temperature metallic state. A phase-diagram 
of the model is presented at the end of this section. In 
Sec. V we study the magnetic properties of the model. Fi- 
nally, we present a magnetic phase-diagram of the model 
which illustrates how magnetic competition helps stabi- 
lizing the incoherent behavior at low temperature. 



II. MODEL AND THEORETICAL 
FRAMEWORK 

A. The Periodic Anderson Model 

The periodic Anderson model describes a set of non- 
dispersive strongly correlated electrons, locally hybridiz- 
ing with a band of conduction electrons. The model 
Hamiltonian is written in the following form: 

H = Hq + Hfp + Hi 



<ij>a 



(1) 



hi=uj: 



The operators Pier (p^) destroy (create) conduction band 
electrons with hopping amplitude Uj and energy Cp. The 
operators fi^ {fta) destroy (create) electrons in the non- 
dispersive orbital with energy e/. The terms propor- 
tional to Vfp describe the hybridization between the two 
species. The interaction term Hi describes the strong 
on-site Coulomb repulsion experienced by /-orbital elec- 
trons. 

The non-interacting lattice Green's function reads: 



-Vfp 



/i - ep - e(k) 



with e(k) the dispersion of the conduction electrons: 
e(k) = 'J2<ij> ^ij- The corresponding inter- 

acting Green's function can be expressed by means of the 
following matrix Dyson equation: 



Go-(k,iw„) ^ = G'o^^(k, iw„) - i;<j(k,ii:j„) (2) 



where 



E/^(k,iw„) 








(3) 



is the self-energy matrix So-- The non- interacting nature 
of the conduction band is reflected in the existence of 
only one non-zero self-energy for the /-electrons. Never- 
theless, it is useful to define an effective self-energy also 
for the conduction electrons as: 



5jp(T(k, iio„ 



-/i- £/ - S/CT(k, iWn) 



(4) 



This function describes the dressing of the p-electrons as 
an effect of both their hybridization with the correlated 
/-electrons and, indirectly, of the Hubbard repulsion on 
the latter. In particular, the appearance of a finite imag- 
inary part in the zero-frequency limit signals the break- 
down of a FL picture for the conduction electrons. 

Since Epo- arises due to both the hybridization and the 
interaction U , it is not expected to vanish in the non- 
interacting limit U = Q. On the other hand, it is easy to 
realize that in this limit the pure hybridization can not 
lead to a finite imaginary part of Epcr at zero frequency, 
and that any breakdown of the FL behavior can descend 
only from correlation effects. 



B. DMFT equations 

The PAM has been studied using a large variety of 
numericaP°"^^ and analytical methods. ^^"^^ To access 
the non-perturbative regime of the PAM, we investigate 
the solution of the model using the DMFT, which has 
been used to solve this model since its early stages. ^^''^'^ 

Within DMFT a lattice model is mapped onto an effec- 
tive single-impurity problem, fixed by a self-consistency 
condition which enforces the equivalence between the two 
models as far as the local physics is concerned. ^^'''^ The 
scheme is equivalent to a local approximation on the self- 
energy, which becomes momentum independent. 

The DMFT equations can be obtained using a quan- 
tum cavity method. The effective action of the single 
/-orbital impurity problem is obtained integrating out 
all lattice degrees of freedom except for a chosen site (la- 
beled conventionally as site i = 0) and keeping only the 
first term in the expansion'^^^'''^ in terms of many-particle 
Green's functions: 

Seff =- f dr f dr' ^ /o+ (r)C7-; (r - r')/o.(r') 
Jo Jo 



+ U / driufotir) - l/2][n;oi(T) - 1/2] 
"'0 



(5) 



The action S'off is expressed in terms of the local 
Weiss Field t/g^^(iwn), describing the quantum fluctu- 
ations at the correlated /-orbital. The Weiss field sat- 
isfies a self-consistence equation which depends on the 
lattice under consideration. In this work we consider a 
Bcthe lattice with semi-elliptical density of states of half- 
bandwidth D (fixing the energy unit of the problem). 
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(a) doped Kondo Insulator (b) doped Mott Insulator 

Figure 1: (Color online) Schematic representation of the 
doped Kondo (a) and doped Mott (b) insulators, inspired by 
Ref. 33. A schematic picture of the band structure is reported 
in the left side of each figure. The hybridized bands are indi- 
cated with thick (blue online) lines, dashed lines are used for 
the bare bands {Vfp = U = 0). In the right sides we sketched 
the corresponding density of states. A hghter color indicates 
the more correlated character of the hybridized bands. 
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\J — e^. In this case the self-consistency 



is particularly simple and reads: 



iuJn) = iw„ -I- ^ - £/ 



V? 
fp 



- /i — ep j-Gpa-{iuJn 



(6) 

where Gpo- is the conduction electron local Green's func- 
tion. The functional form of Qq^ mirrors in the DMFT 
equations the relation between the two orbitals in the 
lattice problem. The fluctuations at the /-orbital are in 
fact composed of two contributions: (a) the on-site quan- 
tum fluctuations and (b) indirect delocalization through 
conduction band proportional to squared hybridization 
amplitude. 

The DMFT solution requires therefore to compute the 
impurity Green's function 



(7) 



where the symbol Os^a indicates the average with re- 
spect to the effective action (5). From the knowledge 
of the impurity Green's function it is straightforward to 
determine the self-energy: 



and finally to evaluate the local Green's function: 



D{e) 



Sp<T(iw„) 



Then a new Weiss field can be computed and the proce- 
dure can be iterated until convergence is achieved. 

The solution of the effective impurity problem, 
i.e. Eq. (7), is the bottleneck of the DMFT algo- 
rithm. In this work we use a combination of nu- 
merical techniques: Hirsch-Fye Quantum Monte Carlo 
(QMC)'"^'*'-''''' and Exact Diagonahzation (ED) methods, 




Figure 2: (Color online) Evolution of the /- (solid line) and p- 
orbital (dashed line) projections of the DOS. Data from QMC 
calculations at T = 0.0125, Vfp — 0.9, ritot ~ 3, analytically 
continued on the real axis using Maximum Entropy Method^^ . 
The figure qualitatively illustrates the Mott metal-insulator 
transition driven by correlation in the PAM. 



both in the full diagonahzation and Lanczos algorithm 
implementations, at zero'^'' and finite temperature.''^ The 
ED method is based on a discretization of the effective 
bath on an adaptive energy grid. In this paper we present 
full ED calculations in which the bath is described by 7 
energy levels and Lanczos calculations with 8 levels. The 
ED calculations have been cross-checked against Density 
Matrix Renormalization Group, which allows to substan- 
tially increase the number of bath levels. 



III. THE HOLE-DOPED MOTT INSULATOR 

The PAM has been largely investigated in proximity 
of the Kondo insulator regime. ''^"''^^ The Kondo insulator 
is a band insulator realized at even integer total filling 
(j^tot = 2). In this regime the system has two hybridized 
bands with a central Kondo peak, corresponding to the 
resonant scattering of the conduction electrons on the 
localized moments and split by an indirect gap Aind (see 
fig. 1(a)). Upon doping the Kondo resonance remains 
pinned at the chemical potential, and the system behaves 
like a heavy-fermion liquid. 

In this work we focus on a different model regime, 
namely the correlated metal obtained by a state with 
odd total occupation (ritot = 1 or 3) and large enough 
interaction. In the case of ritot = 1 or 3, an important 
role is played by the ratio U/A, where A — \ep ~ ef \ is 
the charge-transfer energy, i.e. the separation in energy 
between the two electron orbitals. Two regimes can be 
distinguished:"*^ (a) for A smaller than U the model is 
in the so-called Charge- Transfer (CT) regime, that is ex- 
pected to capture the properties of intermediate to late 
transition-metal oxides. Nevertheless in these systems 
the non-local hybridizations become important and re- 
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5 5 

Figure 3: (Color online) Left panel: renormalization constant 
Z as a function of the doping 5 = |3 — ntot|- Data are from 
Lanczos ED at T = 0, C/ = 2, Vfp = 0.9 and increasing size 
of the effective bath Ns = 8 (diamonds, circles) and 10 (tri- 
angles, squares). Right panel: moment-moment correlation 
— {rrizp ■ rrizf) as a function of the doping 5 in the hole doped 
regime. Data are from full ED calculations (see Appendix A) 
for T — 0.008 and for the same model parameters. 



quire the introduction of other terms in the Hamiltonian 
to be properly described, (b) For A larger than U the 
model is in the Mott-Hubbard (MH) regime, which models 
the properties of early transition-metal oxides and heavy 
fermion systems, usually dominated by local physics. In 
this work we shall focus on this latter model regime and 
study the doping of a Mott insulator. 

In the simplest sketch of this regime, the non- 
correlated band has a lower energy than the correlated 
one (which however is dispersive only because of the hy- 
bridization with the itinerant fermions) . The latter band 
is in turn split by the Mott gap (see fig. 1(b)). Simi- 
larly to the Kondo Insulating regime, a heavy fermion 
state is obtained upon finite doping as soon as the sys- 
tem develops a coherent Kondo resonance signature of 
the insulator-metal transition. 

In the following we shall briefiy review the formation 
of a correlated metallic state by hole doping^^'^^. With- 
out loss of generality, we fix the energy of the correlated 
orbitals at the Fermi level ef = and Cp ~ —1, so that 
A = 1. For U = Vfp = the model describes a system 
with completely filled conduction band Up = 2 and half- 
filled correlated orbitals with nj ~ \. For finite values of 
the hybridization the correlated electrons can move with 
an effective hopping of the order of ieff — V^p/A, corre- 
sponding to the indirect delocalization through the con- 
duction band (see top panel of fig. 2). The hybridization 
modifies the orbitals occupation, pushing a substantial 
amount of the p-electron states to the Fermi level, so that 
Up < 2 and consequently Uf > 1 and the relevant carriers 
are hybrid in nature. However, the /- and p-character of 
the model solution can still be used to indicate the pro- 
jection onto the correlated and non-correlated orbital rc- 
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Figure 4: (Color online) Evolution of the imaginary part of 
the conduction electrons self-energy ImEp(iaj„) for increasing 
temperature. Data are from finite temperature Lanczos ED 
with doping S = 0.1. 



spectively. Upon increasing the interaction strength (see 
central panel of fig. 2) we first observe the formation of 
a correlated metallic state. This is characterized by the 
presence in the density of states (DOS) of a metallic fea- 
ture at the Fermi level, fianked by the two precursors of 
the Hubbard bands. A Mott insulating state is then ob- 
tained further increasing the correlation U. The system 
opens a spectral gap at the Fermi energy (see bottom 
panel of fig. 2) with a width controlled by the correla- 
tion U.^^ To fix ideas, in the remaining part of this work 
we shall set the correlation and the hybridization to, re- 
spectively, U = 2.0 and Vfp = 0.9. This choice of the 
model parameters corresponds to a Mott insulating state 
for ntot = 3. Similar results can be obtained for different 
values of correlation and hybridization. 

The Mott insulator can be destabilized in favor of a 
correlated metallic phase by either adding or removing 
electrons (creating holes). As first noticed in Rcf. 17 the 
two transitions have a different character, ultimately re- 
lated to the different role played by the non interacting 
band in the two cases. Doping with electrons, the ex- 
tra carriers populate essentially the correlated orbitals 
while the p-band remains almost filled and its role is to 
allow the delocalization of correlated electrons. In other 
words, in this regime there are no multi-band effects and 
the hybridization plays a minor role. Therefore the f- 
electrons behave essentially as in a single-band Hubbard 
model with an effective hopping of the order of ioff ■ 

In the hole doped regime the electronic configuration is 
substantially different. The holes are essentially associ- 
ated to the absence of p-electrons, and they tend to bind 
to the local moments of the almost half-filled correlated 
orbitals^^. It is already apparent that this state can not 
be described by a single-band model. This effect is evi- 
dent in the behavior of the inter-orbital moment-moment 
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Figure 5: (Color online) Main panel: imaginary part of 
the conduction electrons self-energy lmSp{i(jj„) for increasing 
value of the hole doping and T = 0.001. Data are from Lanc- 
zos ED calculations. Inset: Comparison of the ImSp{iuj„) 
behavior from different numerical methods for 5 = 0.05. 
The other model parameters are the same as in the main 
panel. The QMC and full ED calculations are performed at 
T = 0.008. DMRG is a T = calculation performed with a 
cluster of Ns = 30 sites and plotted down to the position of 
the lowest energy pole. 



correlation 

reported in fig. 3, which shows how the moment of the 
doped p-holes aligns with the moment of the localized f- 
electron. The doping-driven metalization appears as the 
process of delocalizing a multi-band "Zhang-Rice-like" 
singlet state, formed by an itinerant hole bound to a lo- 
calized spin, similar to that proposed in the framework 
of the high-Tc superconductors. '"^ The low-energy prop- 
erties of this metallic state can not be straightforwardly 
interpreted in terms of a single-band Hubbard model, ^''^^ 
and it leads to remarkable properties. 

A first partial indication of the anomalous nature of 
this state comes from an evaluation of the quasi-particle 
weight Z = [1 ~ dlmT,f{iu})/dLo]~^ ^, which measures 
the degree of metallicity of a system, being zero for a 
Mott insulator and one for a non-interacting metal. The 
results (see fig. 3) show that Z is substantially smaller for 
the hole-doped than for the electron-doped case, already 
signaling that the Zhang-Rice liquid is a poorer metal 
than a standard correlated metal. In the following we 
will show that the difference goes well beyond the quasi- 
particle renormalization. 



IV. THERMAL BREAKDOWN OF THE 
FERMI-LIQUID 

Fermi-liquid theory is the standard paradigm for 
metallic systems and describes correlated Fermi systems 
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Figure 6: a — lmEp(ia;n — 0) as a function of the temper- 
ature for different values of the doping. The data shown are 
from full ED calculations. 



as a collection of non-interacting renormalized quasi- 
particles. DMFT studies of various correlated models 
have shown that even very close to the Mott transi- 
tion the correlated metallic state is typically a Fermi 
liquid with a reduced effective hopping proportional to 
the quasi-particle weight Z. This scale also controls the 
coherence temperature above which the coherent motion 
of the carriers is destroyed by thermal fluctuations. 

In this section we will show that the correlated metallic 
state of the PAM in the weakly hole-doped regime turns 
out to be very fragile with respect to small temperatures. 
More precisely, our system will be a Fermi liquid only 
below an extremely small coherence temperature which, 
for small doping, can be substantially smaller than the 
renormalized Fermi energy controlled by Z. Therefore 
the corresponding metallic state can not be described in 
terms of long-lived quasi-particles but is rather a liquid 
of short-lived singlet-like electronic excitations. 

To substantiate this discussion we study the evolution 
of the imaginary part of the conduction electron self- 
energy IniSp(iw„). The results of our calculations for 
(5 = 0.1 are presented in fig. 4. A Fermi liquid state 
corresponds to a linear behavior of Iml]p(ia;„) at low fre- 
quency, observed only at the lowest investigated tempera- 
ture, T = 0.0005. When we increase T at values of the or- 
der of T = 0.0007, two orders of magnitude smaller than 
the renormalized Fermi energy, the conduction-electron 
self-energy does not vanish in the w — > limit, signaling 
a departure from the Fermi-liquid paradigm. Further in- 
creasing the temperature leads to an enhancement of this 
anomaly. 

In fig. 5 wc follow the evolution of Iml]p(ia;„) for in- 
creasing doping at r = 0.001. For small doping we have a 
clear NFL increase at small frequency which survives up 
to 6 ~ 0.16. For larger doping the system is not strongly 
sensitive to the Mott-Hubbard physics and the standard 
Fermi- liquid behavior is restored around 5 ~ 0.2. 

The violation of the Fermi-liquid paradigm can be 
summarized by the temperature dependence of (t(T) = 
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from Lanczos ED calculations. Lines are guide to the eye. T 



ImEp(iw„ — > 0), reported in fig. 6. This quantity is re- 
lated to the scattering rate of the carriers. In a metallic 
regime ct{T) is expected to vanish at low temperature. 
While for large doping (right panel) a{T) vanishes as 
T — >■ (even if for S = 0.15 some anomaly is observed at 
intermediate temperature), the small-doping data clearly 
confirm the NFL behavior down to very small tempera- 
ture, even if, strictly at T = the vanishing a would be 
recovered. 

Finally, fig. 7 depicts the inverse life-time t^^ — Z^a 
of the doped carriers, where = 1 — ImI]p(ia;i)/7rT'. 
In a Fermi liquid grows as ~ uP' at low temper- 
ature. Our calculations for small doping show a decay 
faster than which strengthens the picture of an inco- 
herent metallic state. Once again, a Fermi-liquid behav- 
ior is established only at extremely low temperatures if 
the doping is small, while the large-doping data recover 
the standard behavior. 

The increasing scattering rate as a function of decreas- 
ing temperature is usually associated to scattering with 
impurities^''. In this spirit, in the following we will in- 
terpret our results as the scattering of the carriers with 
fluctuating local moments. This effect can be understood 
as the results of the competition between the aforemen- 
tioned tendency to form local Zhang-Rice-singlets, driven 
by the hole-doping, and the incoherent nature of the scat- 
terer provided by the /-electron local moments, driven by 
Mott physics. At large doping the increased number of 
available holes of p-type helps the formation of a many- 
body coherent state without breaking the local binding 
with /-moment. This arguments will be substantiated by 
the calculations that we report in the following sections. 



A. The coherence temperature 

The analysis of the self-energy and of the carriers 
lifetime clearly shows the existence of a small doping- 
dependent energy scale associated with the appearance 
of an incoherent metal. We expect this scale to influence 



Figure 8: (Color online) Local spin susceptibility xioc as a 
function of the temperature and increasing value of the hole 
doping. Data are from Lanczos ED (open symbols) and full 
ED (pluses, crosses and stars symbols) calculations. 

also other observables, like the local spin susceptibility: 
Jo 

This quantity describes the response to a local mag- 
netic field and easily discriminates between a Fermi- 
liquid, in which the zero-temperature limit is a constant 
(Pauli susceptibility) , and a paramagnetic Mott insulator 
in which it diverges like 1/T (Curie behavior). 

The results are reported in fig. 8. In the Mott insulat- 
ing state ((5 = 0) the magnetic moments of the localized 
/-electrons essentially behave as free spins, we thus ob- 
tain the typical Curie behavior with a 1/T dependence 
for the spin susceptibility. The slightly hole-doped regime 
does not show the behavior of a standard metal, namely 
Xioc keeps on increasing down to the lowest investigated 
temperature T ~ 10~^D without any sign of saturation. 
The enhancement of the spin susceptibility signals the 
presence of unquenched local moments and can be associ- 
ated to protracted screening effect.'*'' Only for larger dop- 
ing, the susceptibility saturates to large constant value at 
very low temperature. 

The presence of enhanced low-T spin susceptibility 
coexisting with a (bad) metallic behavior substantiates 
the idea that the hole-doped system can be regarded as 
formed by nearly free (incoherent) moments, and an un- 
derlying metallic host formed by the doped holes which 
are prevented from coherently delocalize by local cou- 
pling to /-moments. This interpretation leads us to esti- 
mate the coherence temperature Tcoh from 

Xi;J(T)cxr + r,oh (8) 

We plot the resulting values, obtained with different 
numerical methods, in fig. 9. In the same plot we re- 
port the crossover points estimated from the tempera- 
ture evolution of the imaginary part of the self-energy 
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Figure 9: (Color online) Coherence temperature scale Tcoh 
as extrapolated from the inverse local spin susceptibility 
Xioc (^) ■ T^i^^ extrapolations from different numerical methods 
are found to be in satisfactory agreement. 



Figure 10: (Color online) Phase-diagram of the PAM near the 
Mott insulating state as a function of temperature and hole- 
doping. The diagram is obtained from lmEp{iuj„ — >■ 0). The 
dotted line indicates the crossover temperature scale Tcoh. 



Sp (red crosses). The good agreement of these points 
with the extrapolated data validates the physical inter- 
pretation of the coherence temperature. It is unfortu- 
nately very hard to identify the functional form of the 
coherence temperature due to the smallness of the scale 
involved and the numerical uncertainties. However, the 
data are compatible with an exponential behavior of the 
form Tcoh — -Be""*/^, which has been obtained within the 
1/A^ approximation in the infinite- C/ Kondo limit^^. 

The phase diagram in the doping-temperature plane, 
presented in fig. 10, can help us to summarize the sce- 
nario emerging from our calculations. The diagram re- 
veals the character of the DMFT solution in proximity 
of the Mott insulating state through the behavior of the 
ImEp(icj„ — > 0). Using finite temperature Lanczos ED 
method we investigated a smaller temperature scale with 
respect to that studied in Ref. 18. A large value of the 
lmSp{iujn — ^ 0) testifies a NFL behavior and the results 
clearly show that the highly incoherent state emerges 
from the Mott state and occupies a sizable region of the 
phase diagram. The NFL region is separated from the 
coherent metal by a crossover taking place at Tcoh de- 
fined above, which therefore confirms its meaning as the 
temperature in which the metal loses coherence. 



V. MAGNETIC PROPERTIES 

A. External magnetic field 

We have shown that hole-doping the Mott insulating 
phase of the periodic Anderson model leads to peculiar 
charge carriers, so that the motion of the created p- holes 
occurs through the formation of Zhang-Rice singlets, in 
which the spins of the conduction electrons are anti- 
ferromagnetically correlated with the localized spins. As 
a consequence, we expect that a magnetic field can have 



important and surprising effects on this phase, showing a 
further difference with respect to a standard Fermi liquid. 

In the model regime investigated in this work, the main 
source of magnetism comes from the /-electrons. The 
conduction band is almost completely filled, so that the 
magnetization of the few singly occupied orbitals, favored 
by hole doping, is not expected to contribute significantly 
to the magnetic properties of the system. 

Nevertheless, conduction band electrons can be indi- 
rectly affected by the magnetic polarization of the /- 
orbital moments, through their local binding. To illus- 
trate this point, we show in fig. 11 the evolution of the 
low energy part of Iml]p(ia;„) as a function of a uni- 
form magnetic field B. Apparently the NFL state turns 
into a normal metallic state by the action of an exter- 
nal magnetic field. It is however worth noting that the 
Fermi liquid is recovered for B ~ 0.051?, a huge value 
if compared with experimentally accessible fields. This 
large value is a direct consequence of the large (order 
one) value of V/p, chosen to emphasize the hybridization 
effects and their role in the conduction properties of the 
model. Smaller and more realistic values of this param- 
eter are expected to reduce the critical field by reducing 
the charge fiuctuations at correlated /-orbitals. 

The crossover to a Fermi liquid state driven by exter- 
nal magnetic field is not surprising in light of our anal- 
ysis. Upon increasing the magnetic field a larger and 
larger number of local /-moments are polarized. When 
the moments are aligned, the p-holes can move essentially 
freely in the ferromagnetic background without breaking 
the singlet state with the localized spins. Therefore the 
source of scattering disappears and the metallic state re- 
covers the Fermi-liquid coherence. In other words the po- 
larization of /-orbital local moments allows the conduc- 
tion electron cloud to dynamically screen the correlated 
electrons local moments, dramatically increasing the co- 
herence scale of the system. 
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Figure 11: (Color online) Imaginary part of the majority-spin 
p-electron self-energy for increasing external magnetic field B. 
The data are from QMC solution at T = 0.016 and S = 0.05. 




0.20.40.60.8 1 1.21.4) 

CO 



0.5 



1.5 



(0 



Figure 12: Imaginary part of conduction band electron 
Green's function ImGpcr(icJn) for T — 0.016, 5 = 0.05 and in- 
creasing strength of external magnetic field. Data from QMC 
calculations. 



The coherent motion of the doped carriers with the 
opposite spin of the localized momenta (majority spin) 
should then be balanced by the insulating nature of the 
minority spins carriers. This effect is illustrated in the 
fig. 12. In this figure we show the behavior of both 
spin species conduction electrons Green's function for the 
same strengths of the external magnetic field as used in 
fig. 11. Left panel shows the increasing metalization of 
the majority spin charge carriers, whereas in the right 
panel we show how minority spins are driven towards an 
insulating state by increasing magnetic field. 
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Figure 13: (Color online) Main panel: Staggered magnetiza- 
tion rriAF as a function of temperature and increasing value 
of hole-doping. The data are from full ED calculations. 



hybridization with p-orbital states. In this section we in- 
vestigate the onset of an AFM long-range ordered state 
and its effect on the coherence scale using the extension 
of the DMFT equations to long-range order detailed in 
Appendix B. 

To begin with we report in fig. 13 the staggered magne- 
tization niAF = l/^X]i(^l)'('^/it ^ "/ii) ^ function 
of the temperature for various doping. The transition ap- 
pears to be of second order in the whole doping region. 
The Neel temperature T/v, extracted from a power-law 
fit of the data, is maximum at zero doping and decreases 
by adding holes, as in the single-band Hubbard model. ^'^ 

The onset of an AFM ordering of the local /-moments 
reinstates the Fermi liquid properties in the tiny hole- 
doped regime. This effect is illustrated in the left panel 
of fig. 14, where we present the evolution of the imagi- 
nary part of Tipa-{iuin) from the paramagnetic NFL phase 
to the AFM ordered phase. The large finite intercept 
present in NFL phase is driven to zero in the AFM or- 
dered phase. Nevertheless, the metallic character of the 
solution is preserved across the transition, as illustrated 
in the right panel of the same figure by comparing the 
imaginary parts of the conduction electrons Green's func- 
tions in the two phases. The ordering of the local mo- 
ments in a Neel state, allows the doped charge carriers 
to form coherent electronic waves (with doubled wave- 
vector) and to get delocalizcd. However, as mentioned 
above, the AFM state is only stable in a small window 
of doping and the NFL remains stable for a wide range 
of parameters. 



B. Anti-ferromagnetic ordering 



C. Magnetic stability 



At low temperature, we expect the development of 
anti-ferromagnetic (AFM) correlations a result of super- 
exchange between neighboring /-electrons assisted by the 



The common wisdom about systems of concentrated 
impurities described by the PAM is that long-range mag- 
netic ordering is likely to set in, especially if the metallic 
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Figure 14: (Color online) Conduction electrons self-energy 
luiEpcr{iijJn) (left panel) and Creen's function ImGp(iaj„) 
(right panel). Data from full ED calculations for S = 0.01 
and T = 0.005. 
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Figure 15: (Color online) Intensity plot of ImEp(Lj — >■ 0) as a 
function of external magnetic field B and T at fixed doping 
S = 0.01. For visualization, the data have been normalized to 
max{ImEp(Lj„)} at each (B, T). Dashed lines are drawn to 
better visualize the crossover regions in the phase-diagram. 



state is weakened by correlations as in our case. 

Our results for the AFM state suggest instead a re- 
markable stability of the incoherent metallic state as the 
long-range order is confined to low temperature and small 
doping concentration. In this section we discuss the phys- 
ical origin of this surprising result. 

At small doping near the Mott insulating state neigh- 
boring /-orbital electrons develop AFM correlations as 
a results of super-exchange. These processes are of the 
fourth order in the hybridization with a leading energy 
scale of the orderi^s Jg^ oc W^e^/U ~ V^^/A'^U 

On the other hand it is easy to realize that at large 
doping ferromagnetic correlations are expected because 
of the fact that the doped carriers are locked in singlets 
with the localized /-spins. In an AFM or disordered 
background, the motion of the p-holcs requires to break 
the singlet and it is therefore strongly inhibited, lead- 
ing to the lack of coherence that we discussed at length. 
Moreover it leaves a local moment unscreened, increasing 
the fluctuations in the local magnetization. Conversely, 
a ferromagnetic alignment of the localized spins allows 
for an unperturbed delocalization of the carriers, with 
a mechanism which is closely reminiscent of the double- 
exchange,''^"^^ where the coupling between conduction 
electrons and localized spins is given by the ferromag- 
netic Hund's coupling. 

Therefore, upon increasing the doping the tendency to 
form AFM ordering is contrasted by the increased rel- 
evance of the kinetic energy and eventually it becomes 
more favorable to sacrifice the gain in super-exchange 
energy in order to gain the kinetic energy associated to 
the ferromagnetic background. This leads, most impor- 
tantly, to an intermediate region between the two regimes 
in which the local magnetization is strongly fluctuating. 

From this discussion it is natural to associate the fluc- 
tuations of the local magnetization to the scattering 
mechanism that leads to the poor coherence. To test this 



idea we study the response of the system in the AFM 
ordered metallic phase to the application of a uniform 
magnetic field B which will clearly favor the ferromag- 
netic tendency. 

The results are summarized in the phase diagram of 
fig. 15, determined again using ImEp(iw„ — >■ 0). Details 
of the calculations are given in Appendix B. Doping is 
fixed to S — 0.01, safely into the AFM ordered region in 
the hmit B -> 0. 

The phase-diagram shows that the AFM order sur- 
vives the effects of the external magnetic field up to small 
strengths {B ~ 0.01). In this region the solution keeps 
the coherent metallic character enforced by the long- 
range magnetic ordering. Nevertheless, for larger values 
of the magnetic field the system is driven to an incoherent 
state with finite-temperature NFL behavior, as indicated 
by the increased scattering (light color). In this regime 
the large applied field tends to magnetically polarize the 
AFM ordered /-moments, producing their strong frustra- 
tion and ultimately leading to the formation of an inco- 
herent magnetic background for the motion of the doped 
carriers. Further increasing the strength of the magnetic 
field triggers the formation of a ferromagnetic ordering of 
the /-moments and a fully polarized (coherent) metallic 
state (right dark coloured area). 

The most striking observation is that the present di- 
agram faithfully mirrors the diagram as a function of 
doping, clearly suggesting that the evolution of the con- 
duction properties as a function of doping is associated 
to the transition from the AFM state to the ferromag- 
netic regime and that the poorly coherent metal estab- 
lishes precisely in the intermediate region, dominated by 
the local spin fluctuations which appear as the source of 
the scattering mechanism which opposes to the coherent 
motion of the holes. 
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VI. CONCLUSIONS AND PERSPECTIVES 

In this work we presented a detailed dynamical mean- 
field theory study of the properties of the unconventional 
metallic state obtained by doping with holes the Mott 
insulator in the periodic Anderson model. We discuss 
in details the non-Fcrmi liquid behavior of the system 
and the mechanism that is behind the suppression of the 
coherence scale. 

In this regime the holes have mainly p-character, but 
they tend to bind to the correlated /-electrons to form a 
Zhang- Rice-like singlet state. The formation of this com- 
posite object leads to a highly incoherent metallic state 
which deviates from a standard Fermi-liquid above a co- 
herence temperature which decreases very rapidly upon 
reducing doping, and it is much smaller than the effective 
Fermi energy that one could estimate from the degree of 
correlation of the system, i.e. T^oh ^ ZD. 

We characterize this anomalous behavior by studying 
the scattering properties of the carriers and by computing 
the inverse lifetime and local spin susceptibility, which 
allow us to quantitatively estimate the coherence tem- 
perature characterizing the breakdown of the standard 
Fermi liquid and to describe the onset of an incoherent 
metal with finite lifetime. 

The highly incoherent metal is unstable towards anti- 
ferromagnetic ordering only at very small doping, while 
at large doping ferromagnetic correlations develop and 
favour a regular metallic behavior supported by a mech- 
anism which reminds the double-exchange physics. The 
intermediate region, where the motion of the holes is not 
coherent, is therefore dominated by large fiuctuations of 
the /-spins, which provide the scattering channel respon- 
sible of the finite lifetime of the carriers. 

The relation between magnetic fluctuations and the 
breakdown of the standard FL scenario is emphasized by 
observing that an external uniform magnetic field, which 
obviously destroys AFM ordering favoring a ferromag- 
netic alignment, mirrors the effect of doping and leads 
again to a wide region of high incoherence between the 
two magnetically ordered states. 

We emphasize that the path to poor coherence dis- 
cussed in this paper only depends on two general fea- 
tures of strongly correlated materials, namely the Mott 
physics which leads to the localization of carriers and 
multi-orbital physics necessary to the local singlet for- 
mation. In this light, we expect that the mechanism 
outlined here can be a rather general source of violation 
of Fermi liquid paradigm and incoherent behavior, and it 
can be relevant for example to heavy fermions, but also, 
with some important differences related to the d-wave 
symmetry of the Zhang-Rice singlets, to the cuprate su- 
perconductors. 

Finally, a natural question to address is to what extent 
our findings can be considered the local portrait of the 
presence of a quantum critical point, hidden by the ab- 
sence of spatial fluctuations. Indeed, the existence of a 
quantum critical point in the PAM, although in a differ- 



ent model regime, has already been pointed in Ref. 52, 
using cluster extension of the DMFT. The development of 
our work along this direction, in order to clarify the fate 
of the small coherence scale in presence of short-range 
spatial fluctuations, is left for future research. 
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Appendix A: Two-orbital effective impurity model 

The calculation of physical quantities internal to the lo- 
cal "p/-dimer", such as the moment- moment correlation 
function {nizp ■ rrizf) can be performed within single-site 
DMFT using an alternative formulation of the effective 
impurity problem in which the local p-orbital is not in- 
tegrated out in the construction of the effective action. 
Thus the original lattice system is reduced to the prob- 
lem of a single dimer embedded in an electronic bath. 
The corresponding effective action has a 2 x 2 matrix 
structure in the orbital space and reads: 

^cff ^- f dr f dr' J2 i'Lir)Go\r - T')V^o.(r') 
Jo Jo 

+ U / dT[nfot{r) - l/2]Km(r) - 1/2] 
Jo 

The Weiss Field Go^i'i'^n) describes the local quantum 
fluctuations at the tagged dimer. The Bethe lattice self- 
consistency becomes 

GnHiu^n) = ( + ^ - «P - GpiiUJn) -Vfp \ 

The DMFT algorithm for the two-orbital representa- 
tion proceeds as in the standard casse. The effective two- 
orbital impurity problem is solved to determine the im- 
purity Green's functions: 

with a = p, /. Next, the conduction electron self-energy 
Ep can be determined using the Dyson equation and used 
to evaluate the local Green's function Gp which is neces- 
sary to update the local Weiss fleld. The whole algorithm 
is iterated until convergence is reached. 
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Appendix B: Long range order 



and thus: 



The DMFT equations can be extended to describe 
phases with long range magnetic ordering^®. Here we 
derive the equations for the anti-ferromagnetic order in 
the two-orbital effective problem, considering also the ef- 
fect of a uniform magnetic field. Similar equations can 
be derived for the single-orbital effective model. 

On a bipartite lattice crystal as our Bethe lattice, we 
can define two sub-lattices A and B, such that nearest- 
neighbor hopping always connects one ^-site with a B- 
site. Then wc can introduce a four-component spinor 
with orbital and sub-lattice indices so that the bare lat- 
tice propagator takes the form: 



G. 



Okcr 



/ aA -e(k) -Vfp 
-e(k) ae 

-Vfp iujn - ef + HA 

\ ^^/p iuj„ 









with as = ioJn — £p + fJ-s and s = A,B. The corre- 
sponding Green's functions are obtained via the Dyson 
equation with the diagonal self-energy matrix with com- 
ponents {0, 0, E^CT, Sscr}. The p-electrons local Green's 
functions, required to close the DMFT equations, now 
read: 



GpAa{Ki^n) = X! 



e(k)2 



C,Aa — Cb-ct — Co- 



(Bl) 



wc can eliminate one of the two sublattices and recover 
a 2 X 2 formalism with a Weiss field given by 



Gp-cr{i^n) 



-Vfp 



The local conduction electron Green's function 
GpaiiuJn) can be expressed in terms of the following 
Hilbert transform: 



GpaiiuJn) = C- 



de- 



which closes the set of DMFT equations. 

In presence of a uniform magnetic field B in the 
ordered phase of the system, the symmetry relation 
Eq. (Bl) between the two sublattices does not hold. 
Therefore the DMFT solution requires to explicitly con- 
sider the two sublattices and the self-consistency equa- 
tions for the four components of the Weiss field Goas(i^n) 
(s = A, i? and s = B,A) read 



where: 



V? 

fp 



V? 

fp 



^Gpasiil^n) 



7s 



In the case of anti-fcrromagnctic ordering it is not nec- 
essary to take explicitly into account both sublattices. 
Observing that: 



where the coupling to the magnetic field B has been in- 
cluded in a redefinition of the chemical potential Jl^^ = 
MsCT + o'B/2. This means that at each iteration we need 
to solve two impurity models, one for each sub-lattice 
and that the solution of one sub-lattice will determine 
the Weiss field for the other. 
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